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Quantum interferometric optical lithography: 
towards arbitrary two-dimensional patterns 
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As demonstrated by Boto et al. [Phys. Rev. Lett. 85, 2733 (2000)], quantum lithography offers 
an increase in resolution below the diffraction limit. Here, we generalize this procedure in order to 
create patterns in one and two dimensions. This renders quantum lithography a potentially useful 
tool in nanotechnology. 

PACS numbers: 42.50.Hz, 42.25.Hz, 42.65.-k, 85.40.Hp 

Optical lithography is a widely used printing method. 
In this process light is used to etch a substrate. The ex- 
posed or unexposed areas on the substrate then define 
the pattern. In particular, the micro-chip industry uses 
lithography to produce smaller and smaller processors. 
However, classical optical lithography can only achieve a 
resolution comparable to the wavelength of the light used 
H||. It therefore minimizes the scale of the patterns. 
To create smaller patterns we need to venture beyond 
this classical boundary [|J. In Ref. |jj we introduced a 
procedure called quantum lithography that offers an in- 
crease in resolution beyond the diffraction limit. This 
process allows us to write closely spaced lines in one di- 
mension. However, for practical purposes (e.g., optical 
surface etching) we need to create more complicated pat- 
terns in both one and two dimensions. Here, we study 
how quantum lithography can be extended to create these 
patterns. 

This paper is organized as follows: first, for complete- 
ness, we present a derivation of the Rayleigh diffraction 
limit. Then, in Sec. [n]we reiterate the method introduced 
in Ref. ||. Then, in Sec. HI we give a generalized ver- 
sion of the states used in this procedure. We show how we 
can tailor arbitrary one-dimensional patterns with these 
states. In Sec. IV we show how four-mode entangled 
states lead to patterns in two dimensions. Sec. [v] ad- 
dresses the physical implementation of quantum lithog- 
raphy. 



Suppose two plane waves characterised by k\ and ki 
hit a surface under an angle from the normal vector. 
The wave vectors are given by 

fci = fc(cos 9, sin 0) and &2 = k(cos9, — s'm8) , (1) 

where we used |fci| = Ifel = k. The wave number k 
is related to the wavelength of the light according to 
k = 2n/\. 

In order to find the interference pattern in the inten- 
sity J, we sum the two plane waves at position r at the 
amplitude level: 



1(f) oc 



4cos^ 



(2) 



When we calculate the inner product (fci — A^) ■ r/2 from 
Eq. (fj]) we obtain the expression 



I(x) oc cos 2 (fcxsin( 



(3) 



for the intensity along the substrate in direction x. 

The Rayleigh criterion states that the minimal resolv- 
able feature size Ax corresponds to the distance between 
an intensity maximum and an adjacent minimum. From 
Eq. (0) we obtain 



fcAx sin 



(4) 



This means that the maximum resolution is given by 



I. CLASSICAL RESOLUTION LIMIT 

When we talk about optical resolution, we can mean 
two things: it may denote the minimum distance be- 
tween two nearby points which can still be resolved with 
microscopy. Or it can denote the minimum distance sep- 
arating two points which are printed using lithography. 
In the limit of geometric optics these resolutions would be 
identical. In this section we derive the classical resolution 
limit for interferometric lithography using the so-called 
Rayleigh criterion ||. 



Ax 



A 



2fcsin6> 2 (x sin6 4s ' m0 ' 



(5) 



where A is the wavelength of the light. The maximum res- 
olution is therefore proportional to the wavelength and 
inversely proportional to the sine of the angle between 
the incoming plane waves and the normal. The resolu- 
tion is thus maximal (Ax is minimal) when sin 9 = 1 , or 
9 = 7r/2. This is the grazing limit. The classical diffrac- 
tion limit is therefore Ax = A/4. Note that this deriva- 
tion does not use the approximation sin 9 ~ 9, which is 
common when considering diffraction phenomena. 
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II. INTRODUCTION TO QUANTUM 
LITHOGRAPHY 



In this section we briefly reiterate our method of Rcf. 
fl. It exploits the physical properties of multi-photon 
absorption of a substrate. Suppose we have two inter- 
secting light beams a and b. We place a substrate sen- 
sitive to iV-photon absorption at the position where the 
two beams meet, such that the interference pattern is 
recorded. For simplicity, we consider the grazing limit 
in which the angle 9 off axis for the two beams is ir/2 
(see Fig. |n|). Classically, the interference pattern on the 
substrate has a resolution of the order of A/4, where A 
is the wavelength of the light. However, by using entan- 
gled photon- number states (i.e., inherently non-classical 
states) we can increase the resolution well into the sub- 
wavelength regime. 

How does quantum lithography work? Let the two 
counter-propagating light beams a and b be in the com- 
bined entangled state of N photons 

\i>N)ab=(\N,0) ab + e iN ^\0,N) ab )/V2, (6) 

where tp = kx/2, with k = 2ir/A. We define the mode op- 
erator e = (a + b)/y/2 and its adjoint e' = (a* + w)/^/2. 
The deposition rate A on the substrate is then given by 

A N = (V>jv|<MV>w) with 5 N = — — — , (7) 

i.e., we look at the higher moments of the electric field 
operator j7] |{|. The deposition rate A is measured in 
units of intensity. Leaving the substrate exposed for a 
time t to the light source will result in an exposure pat- 
tern P(<p) = A {ft. After a straightforward calculation 
we see that 

A N cx (1 + cosNtp) . (8) 

We interpret this as follows. A path-differential phase- 
shift <p in light beam b results in a displacement x of the 
interference pattern on the substrate. Using two classi- 
cal waves, a phase-shift of 2n will return the pattern to 
its original position. However, according to Eq. (^), one 
cycle is completed after a shift oi2ir/N. This means that 
a shift of 2ir will displace the pattern N times. In other 
words, we have N times more maxima in the interference 
pattern. These need to be closely spaced, yielding an ef- 
fective Rayleigh resolution of Ax = A/ AN, a factor of N 
below the classical interferometric result of Ax = A/4 [0. 



substrate 

Fig. 1. Two light beams a and b cross each other at the sur- 
face of a photosensitive substrate. The angle between them 
is 29 and they have a relative phase difference tp = kx/2. We 
consider the limit case of 8 — > tv/2. 

Physically, we can interpret this result as follows: in- 
stead of having a state of N single photons, Eq. (||) de- 
scribes an A-photon state. Since the momentum of this 
state is N times as large as the momentum for a single 
photon, the corresponding DeBroglie wavelength is N 
times smaller. The interference of this A-photon state 
with itself on a substrate thus gives a periodic pattern 
with a characteristic resolution dimension of Ax = A/4 N. 



III. GENERAL PATTERNS IN ID 

So far, we have described a method to print a simple 
pattern of evenly spaced lines of sub-wavelength resolu- 
tion. However, for any practical application we need the 
ability to produce more complicated patterns. To this 
end, we introduce the state 

\ipNm)ab = U mv \N-m,m) ab 

+e i ( N - m tee ie ™\m,N-m) ab yV2. (9) 

This is a generalized version of Eq. (^|). In particular, Eq. 
(^) reduces to Eq. (^) when m — and 6 m = 0. Note that 
we included a relative phase e l9m , which will turn out to 
be crucial in the creation of arbitrary one-dimensional 
patterns. 

We can calculate the deposition rate again according to 
the procedure in Sec. 0. As we shall see later, in general, 
we can have superpositions of the states given by Eq. (||). 
We therefore have to take into account the possibility of 
different values of m, yielding a quantity 

= (V>iVm|WiVm') ■ (10) 

Note that this deposition rate depends not only on the 
parameter tp, but also on the relative phases 9 m and 9 m >. 
The deposition rate then becomes 



A Nm ' nc 



N\ ( N 
m/\m' 



„—i(N—m—m')(p—i6 m , — i(m — 7n)tpA(9 m i — m ) 

O O I C> O j 



(11) 
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where ( ) means N\/(N — m)\m\. Obviously, 

{i>Nm\h\^N'mf) = when I {N, N'}. For m = m', 
the deposition rate takes on the form 

&Nm ex Q {1 + cos[(A - 2m)^ + 8 m }} , (12) 

which, in the case of m = and 9 m = 0, coincides with 
Eq. When 8 m is suitably chosen, we see that we 

also have access to deposition rates (1 — cos Nip) and 
(1 ± sin Nip). Apart from this extra phase freedom, 
Eq. ([12]) does not look like an improvement over Eq. (||), 
since N — 2m < N, which means that the resolution de- 
creases. However, we will show later how these states can 
be used to produce non-trivial patterns. 

First, we look at a few special cases of 8 m and m i. 



When we write Aj^™ 



A^'(0,0)occos( 

A#™'(0,7T) cxcos 

A^'(7r,0)ocsin( 
A^>,7T)cxsin( 
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) we have 
(N -2m' 



(13a) 
(13b) 
(13c) 
(13d) 



These relations give the dependence of the matrix ele- 
ments A^™ on 6 m and 8 m i in a more intuitive way than 
Eq. (|TT|) does. Finally, when 6 m = B m * = 8 we obtain 



Ajvm « cos 



(N - 2m)if 



(N - 2m')ip - 6 



(14) 



So far we have only considered generalized deposition 
rates given by Eq. with special values of their pa- 
rameters. We will now turn our attention to the problem 
of creating more arbitrary patterns. 

Note that there are two main, though fundamentally 
different, ways we can superpose the states given by Eq. 
(0) . We can superpose states with different photon num- 
bers n and a fixed distribution m over the two modes: 



N 



|*m) = ^2 a n \ij) nm ) 



(15) 



with a n complex coefficients. This is a superposition of 
states with different total photon number in each branch. 
Alternatively, we can superpose states with a fixed pho- 
ton number TV, but with different distributions m: 



LJV/2J 

|*jv) = ^ a m \ip Nm ) 

rn—0 



(16) 



where [N/2\ denotes the largest integer I with I < N/2 
and a m again the complex coefficients. Every branch in 
this superposition is an A-photon state. 

These two different superpositions can be used to tai- 
lor patterns which are more complicated than just closely 
spaced lines. We will now study these two different meth- 
ods. 



A. The Pseudo-Fourier Method 

The first method, corresponding to the superposition 
given by Eq. ( |l5|) , we will call the pseudo-Fourier method 
(this choice of name will become clear shortly). When we 
calculate the deposition rate A m according to the state 
I "I'm) we immediately see that branches with different 
photon numbers n and n' do not exhibit interference: 

N N 

A m = ^ \ a n\ 2 ''{i>nm\5n\i>nm) = ^2 \ a n\ 2 ^nm • (17) 
n=0 n=0 

Using Eq. jl^ ) the exposure pattern P{f>) — A m t be- 
comes 



N 



P(f>) =t^2c n {l + cos[(n-2m)(p + e n }} 



(18) 



where t is the exposure time and the c n are real and 
positive. Since m < n and m is fixed, we have m = 0. 
We will now prove that this is a Fourier series up to a 
constant. 

A general Fourier expansion of p{f>) can be written as 



N 



P((f) = ^2 ( a n C0S nL P + ^« Sm nl P) 



(19) 



n=0 



Writing Eq. (|1S|) as 



N 



N 



P(f) — t ^2 c n + t *^2 C ™ cos ( nl P + On) 



(20) 



n=0 



where t X^=o c « ^ s a constant. If we ignore this constant 
(its contribution to the deposition rate will give a general 
uniform background exposure of the substrate, since it is 
independent of if) we see that we need 

c„ cos(nip + 9 n ) — a n cosnif + b n sinnif (21) 

with c„ positive, 8 n € [0, 2tt) and a n , b n real. Expand- 
ing the left-hand side and equating terms in cosnif and 
sin nip we find 



c„ cos f 



and b n — c n sin 9 n 



(22) 
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This is essentially a co-ordinate change from Cartesian 
to polar co-ordinates. Thus, Eq. (18) is equivalent to a 
Fourier series up to an additive constant. Since in the 
limit of N — > oo a Fourier series can converge to any 
well-behaved pattern P(<p), this procedure allows us to 
approximate arbitrary patterns in one dimension (up to 
a constant). It is now clear why we call this procedure 
the pseudo-Fourier method. 

However, there is a drawback with this procedure. The 
deposition rate A is a positive definite quantity, which 
means that once the substrate is exposed at a particular 
Fourier component, there is no way this can be undone. 
Technically, Eq. (jl^) can be written as 

N 

P(ip) = Q ■ t + t (a n cos ntp + b n sin rup) , (23) 

n=0 

where Q is the uniform background 'penalty exposure 
rate' Q — X^=o c « we mentioned earlier. The second 
term on the right-hand side is a true Fourier series. Thus 
in the pseudo-Fourier method there is always a minimum 
exposure of the substrate. Ultimately, this penalty can be 
traced to the absence of interference between the terms 
with different photon number in Eq. (^). Next, we will 
investigate whether our second method of tailoring pat- 
terns can remove this penalty exposure. 



B. The Superposition Method 

We will now study our second method of tailoring pat- 
terns, which we call the 'superposition method' (lacking a 
better name). Here we keep the total number of photons 
N constant, and change how the photons are distributed 
between the two beams in each branch [see Eq. ([nil)]. 
A distinct advantage of this method is that it does ex- 
hibit interference between the different branches in the 
superposition, which eliminates the uniform background 
penalty exposure. 

Take for instance a superposition of two distinct terms 



\^n) = a m \ip Nm ) + a m ,\4i Nm i) , 



(24) 



with |a m | 2 + |aw| 2 = 1 and \ip nm ) given by Eq. (]T|). Af- 
ter some algebraic manipulation the deposition rate can 
be written as 



Aw tx la 



N 



{l + cos[(7V- 2m)ip + 9 m ]} 
+ \u m ,\ 2 ( N , ) {1 + cos[(N - 2ro> + m ']} 



+8r. 




x cos - [(N — 2m)tp + 6 
x cos - [(N - 2m')ip + i 



where the deposition rate A is now a function of a m and 
a m ' , where we have chosen the real numbers r™ and £™ 
to satisfy u* m u m ' = r™ exp(i£™ ). For the special values 



N = 20, m = 9, 



= 5 and B m = 9 m i = we obtain the 



pattern shown in Fig. [Tv| . Clearly, there is no uniform 
background penalty exposure here. 



with m = 5 and m! 




Fig. 2: The deposition rate due to a superposition of two 
states containing 20 photons with distributions m = 9 and 
m! = 5 (6m ~ 6m' = 0). The deposition rate at tp — n/2 
and tp = 37r/2 is zero, which means that there is no general 
uniform background exposure. 

For more than two branches in the superposition this 
becomes a complicated function, which is not nearly as 
well understood as a Fourier series. The general expres- 
sion for the deposition rate can be written as 



Ajv oc 



L7V/2J LJV/2J 

EE-i 

m— m'— 
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x cos - [(N 
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where we have chosen 



2m)tp + 6 
2m')tp + 1 

and 



m>] , (26) 
real to satisfy 



(25) 



a* m a m , = C exp(iC ). Note that f™ = 0. 

If we want to tailor a pattern F(<p), it might be the 
case that this type of superposition will also converge to 
the required pattern. We will now compare the superpo- 
sition method with the Fourier method. 



C. Comparing the two methods 

So far, we discussed two methods of creating non- 
trivial patterns in one dimension. The Fourier method 
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is simple but yields a uniform background penalty expo- 
sure. The superposition method is far more complicated, 
but seems to get around the background exposure. Be- 
fore we make a comparison between the two methods we 
will discuss the creation of 'arbitrary' patterns. 

It is well known that any sufficiently well-behaved pe- 
riodic function can be written as an infinite Fourier series 
(we ignore such subtleties which arise when two functions 
differ only at a finite number of points, etc.). However, 
when we create patterns with the pseudo-Fourier lithog- 
raphy method we do not have access to every component 
of the Fourier expansion, since this would involve an in- 
finite number of photons (n — > oo). This means that we 
can only employ truncated Fourier series, and these can 
merely approximate arbitrary patterns. 

The Fourier expansion has the nice property that when 
a series is truncated at N, the remaining terms still give 
the best Fourier expansion of the function up to N. In 
other words, the coefficients of a truncated Fourier series 
are equal to the first N coefficients of a full Fourier series. 
If the full Fourier series is denoted by F and the trun- 
cated series by Fn, we can define the normed-distance 
quantity D^: 



not the same) which can be written as 



D 



N 



2tt 



\F{ip)-F N {ip)\ 2 dip 



(27) 



which can be interpreted as a distance between F and 
Fn- If quantum lithography yields a pattern pn(<p) = 
A^t, we can introduce the following definition: quan- 
tum lithography can approximate arbitrary patterns if 



2 IT 



\F(ip)-P N (ip)\ 2 dip<eD N 



(28) 



with e some real, positive definite proportionality con- 
stant. This definition gives the concept of approximating 
patterns a solid basis. 

We compare the Fourier and the superposition method 
for one special case. We choose the test function 



h\i - \<f<\ 
otherwise . 



(29) 



With up to ten photons, we ask how well the Fourier and 
the superposition method approximate this pattern. 

In the case of the Fourier method the solution is im- 
mediate. The Fourier expansion of the 'trench' function 
given by Eq. (E9T) is well known: 



q=0 



(-1)" 

2q+ 1 



cos[(2g + l)ip] 



(30) 



Using up to n = 10 photons we include terms up to q = 4, 
since 2q + 1 < 10. The Fourier method thus yields a pat- 
tern P{<f) (the two patterns P{f) and F(ip) are generally 



P(<p) 
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E 2q~ 



Cat 



(l + cos[(2g+lV + 7TK 9 ]) , (31) 



where c q is a constant depending on the proportionality 
constant of A 2q +i, the rate of production of \ip n m) an d 
the coupling between the light field and the substrate. 
The term n q is defined to accommodate for the minus 



signs in Eq. (p0|): it is zero when q is even and one when 
q is odd. Note the uniform background penalty exposure 
rate Yln=Q c ql (2?+l)< The result of this method is shown 
in Fig. nVL 

Alternatively, the superposition method employs a 
state 



I* 



N 



LAT/2J 
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m=0 



Ni. 



(32) 



The procedure of finding the best fit with the test func- 
tion is more complicated. We have to minimize the abso- 
lute difference between the deposition rate An (a) times 
the exposure time t and the test function F(<p). We have 
chosen a = (ag, . . . , a n / 2 )- Mathematically, we have to 
evaluate the a and t which minimize d^' 



d N = 



2tv 



\F(<p)- A N (a)t\ 2 d<p , 



with 



A N (a) = (^n\S n \^n) 



(33) 



We have to fit both t and a. Using a genetic optimaliza- 
tion algorithm jl0| (with h = 1, a normalized height of 
the test function) we found that the deposition rate is ac- 
tually very close to zero in the interval 7r/2 < if < 37r/2, 
unlike the pseudo-Fourier method, where we have to pay 
a uniform background penalty. This result implies that 
in this case a superposition of different photon distribu- 
tions m, given a fixed total number of photons N, works 
better than a superposition of different photon number 
states (see Fig. IV). In particular, the fixed photon num- 



ber method allows for the substrate to remain virtually 
unexposed in certain areas. 



P(f) = A 10 (f)t 




3^72 2tt 
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Fig. 3: The deposition rate on the substrate resulting from a 
superposition of states with N — 10 and different m (black 
curve), and also resulting from a superposition of states with 
different n with m — (grey curve). The coefficients of the 
superposition that yield the black curve are optimized us- 
ing a genetic algorithm while the grey curve is a trun- 
cated pseudo-Fourier series. Notice the 'penalty' (displace- 
ment from zero) of the deposition rate for the pseudo-Fourier 
series between 7r/2 and 3n/2. 

We stress that this is merely a comparison for a spe- 
cific example, namely that of the trench target function 
F(ip). We conjecture that the superposition method can 
approximate other arbitrary patterns equally well, but 
we have not yet found a proof. Besides the ability to 
fit an arbitrary pattern, another criterion of comparison 
between the pseudo-Fourier method and the superposi- 
tion method, is the time needed to create the TV-photon 
entangled states. 

Until now, we have only considered sub-wavelength 
resolution in one direction, namely parallel to the direc- 
tion of the beams. However, for practical applications we 
would like sub-wavelength resolution in both directions 
on the substrate. This is the subject of the next section. 



Fig. 4: Four light beams a, 6, c and d cross each other at 
the surface of a photosensitive substrate. The angles between 
a and b and c and d are again taken in the grazing limit of 
9 — tt/2. The relative phase difference between a and b is tp 
and the relative phase difference between c and d is X- 

We can now do the same for the y-direction, employ- 
ing two counter-propagating beams (c and d) in the y- 
direction (see Fig. ). The same conditions apply: we 
consider the limit where the spatial angle 9 off axis ap- 
proaches 7r/2, thus grazing along the substrate's surface. 

Consider the region where the four beams a, b, c and 
d overlap. For real lithography we have to take into ac- 
count the mode shapes, but when we confine ourselves to 
an area with side lengths A (where A is the wavelength 
of the used light) this problem does not arise. 

The class of states on modes a to d that we consider 
here are of the form 



e imip \N -m,m; 0,0) 

+ e* N - m ^e Km \m, N-m;Q, 0) 
+ e ik9 \0,0;N-k,k) 

+ e tiN - k)x e^ k \0,0;k,N - k) 



(34) 



IV. GENERAL PATTERNS IN 2D 

In this section we study how to create two-dimensional 
patterns on a suitable substrate using the quantum 
lithography techniques developed in the previous sec- 
tions. As we have seen, the phase shift tp, in the setup 
given by Fig. EL acts as a parametrization for the deposi- 
tion rate in one dimension. Let's call this the x-direction. 




substrate 



where Cm and C~fc are two relative phases. This is by no 
means the only class of states, but we will restrict our 
discussion to this one for now. Observe that this is a su- 
perposition on the amplitude level, which allows destruc- 
tive interference in the deposition rate in order to create 
dark spots on the substrate. Alternatively, we could have 
used the one-dimensional method [with states given by 
Eq. (^)] in the x- and ^-direction, but this cannot give 
interference effects between the modes a, b and c, d. 

The phase-shifts tp and x m the light beams b and 
d (see Fig. ) result in respective displacements x and y 
of the interference pattern on the substrate. A phase- 
shift of 2-7T in a given direction will displace the pat- 
tern, say, N times. This means that the maxima are 
closer together, yielding an effective resolution equal to 
Ax = Ay = \/AN. This happens in both the x- and the 
y-direction. 

We proceed again as in Sec. [0] by evaluating the N th or- 
der moment <5jv of the electric field operator [see Eq. M)] 
This gives the deposition rate A^™ fc fe 
[with \ii) k Nm ) given by Eq. (|34 
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For the special choice of to' = to and k! — k we have 

2 



4m M 1 + ^[(^ ~ 2m )^ + Cm]) + J (1 + CO S [(N - 2k) X + Cfc]) 



+4 



(to)(P c °4^-^ + ^-^ 

x cos - [(N - 2m)<p - Cm] cos - [(N - 2k) x - Cfc] ■ 



(35) 



(36) 



We can again generalize this method and use superpo- 
sitions of the states given in Eq. ([mJ). Note that there 
are now three numbers N, to and k which can be var- 
ied. Furthermore, as we have seen in the one-dimensional 
case, superpositions of different n do not give interference 
terms in the deposition rate. 

Suppose we want to approximate a pattern F(ip, x), 
with {ip, x} € [0, 2n]. This pattern can always be writ- 
ten in a Fourier expansion: 

oo 

F( ( P>X) = E a pq cos p<f cos qx + b pq cos pip sin qx*. 

p,q=0 

c pq sin pep cos qx + d pq sin pep sin qx ■ (37) 



with a pq , b pq , 



z pq and dpq real. In the previous section 



we showed that quantum lithography could approximate 
the Fourier series of a one-dimensional pattern up to a 
constant displacement. This relied on absence of interfer- 
ence between the terms with different photon numbers. 
The question is now whether we can do the same for 
patterns in two dimensions. Or alternatively, can gen- 
approximate the 



eral superpositions of the state \ip% ri 



pattern F(ip, x) 



From Eq. (35) it is not obvious that we can obtain the 
four trigonometric terms given by the Fourier expansion 
ofEq. M): 



patterns like the one given in Fig. can be inferred from 
Eq. (|35|). In order to find the best fit to an arbitrary 
pattern one has to use a minimization procedure. 

For example, we calculate the total deposition rate due 
to the quantum state | ^ at) , where 



LJV/2J [N/2] 

1*^) = E a mkH N m) 

m=0 fc=0 



(39) 



Here, ot m k are complex coefficients. We now proceed by 
choosing a particular intensity pattern F(ip, x) and opti- 
mising the coefficients a mk for a chosen number of pho- 
tons. The deposition rate due to the state \^>n) is now 



LJV/2J \_N/2\ 

A N(d) = E E ^ nk a mlk ,A N N Zk , (40) 

m,m'—0 k,k'—0 



with a — (o!o,o> Q!o,i • • ■ 7 a N/2,N/2)- We again have to 
evaluate the a and t which minimize 



A cx cos pip cosqx , 

A cx cos pip sinqx , 

A cx sin pip cos qx , 

A cx sin pip sinqx . 



(38a) 
(38b) 
(38c) 
(38d) 



\F(ip,x)-A N (d)t\ 2 dipd X . (41) 



We can therefore not claim that two-dimensional quan- 
tum lithography can approximate arbitrary patterns in 
the sense of one-dimensional lithography. Only simple 



The values of a and t can again be found using a genetic 
algorithm. 
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expression. In that case we only need to determine this 
elementary set of patterns. 

Furthermore, we have to study whether the uniform 
background penalty exposure really presents a practical 
problem. One might argue that a sufficient difference 
between the maximum deposition rate and the uniform 
background penalty exposure is enough to accommodate 
lithography. This depends on the details of the sub- 
strate's reaction to the electro magnetic field. 

Before quantum lithography can be physically imple- 
mented and used in the production of nano circuits, these 
issues have to be addressed satisfactorily. 

VI. CONCLUSIONS 



Fig. 5: A simulation of a two-dimensional intensity pattern 
on an area A 2 , where A denotes the wavelength of the used 
light. Here we modelled a square area with sharp edges. The 
pattern was generated by a Fourier series of up to ten photons 
(see also Fig. 3 for the one-dimensional case). 



V. PHYSICAL IMPLEMENTATION 

With current experimental capabilities, the physical 
implementation of quantum lithography is very challeng- 
ing. In particular, there are two major issues to be dealt 
with before quantum lithography can become a mature 
technology. First of all, we not only need the ability to 
create the entangled photon states given by Eqs. (^) and 
but we should also be able to create coherent su- 
perpositions of these states. One possibility might be to 
use optical components like parametric down-converters. 
Contrary to the results of Ref . jll| , we are not concerned 
with the usually large vacuum contribution of these pro- 
cesses, since the vacuum will not contribute to the spatial 
profile of the deposition [see Eqs. @ and (fjp]. 

Secondly, we need substrates which are sensitive to the 
higher moments of the electric field operator. When we 
want to use the pseudo-Fourier method, up to N photons 
for quantum lithography in one dimension, the substrate 
needs to be reasonably sensitive to all the higher mo- 
ments up to N, the maximum photon number. Alterna- 
tively, we can use the superposition method for N pho- 
tons when the substrate is sensitive to predominantly one 
higher moment corresponding to N photons. Generally, 
the method of lithography determines the requirements 
of the substrate. 

There are also some considerations about the approxi- 
mation of patterns. For example, we might not need arbi- 
trary patterns. It might be the case that it is sufficient to 
have a set of patterns which can then be used to generate 
any desired circuit. This is analogous to having a univer- 
sal set of logical gates, permitting any conceivable logical 



In this paper we have generalized the theory of quan- 
tum lithography as first outlined in Ref. || . In particular, 
we have shown how we can create arbitrary patterns in 
one dimension, albeit with a uniform background penalty 
exposure. We can also create some patterns in two di- 
mensions, but we have no proof that this method can be 
extended to give arbitrary patterns. 

For lithography in one dimension we distinguish two 
methods: the pseudo-Fourier method' and the superpo- 
sition method. The pseudo-Fourier method is concep- 
tually easier since it depends on Fourier analysis, but it 
also involves a finite amount of unwanted exposure of the 
substrate. More specifically, the deposition rate equals 
the pattern in its Fourier basis plus a term yielding un- 
wanted background exposure. The superposition method 
gets around this problem and seems to give better results, 
but lacks the intuitive clarity of the Fourier method. Fur- 
thermore, we do not have a proof that t his me thod can 
approximate arbitrary patterns (see Sec. Ill C for a dis- 
cussion on this approximation). 

Quantum lithography in two dimensions is more in- 
volved. Starting with a superposition of states, given 
by Eq. (jU), we found that we can indeed create two- 
dimensional patterns with sub-wavelength resolution, 
but we do not have a proof that we can create arbitrary 
patterns. Nevertheless, we might be able to create a cer- 
tain set of elementary basis patterns. 

There are several issues to be addressed in the future. 
First, we need to study the specific restrictions on the 
substrate and how we can physically realize them. Sec- 
ondly, we need to create the various entangled states in- 
volved in the quantum lithography protocol. 

Finally, G.S. Agarwal and R. Boyd have called to our 
attention that quantum lithography works also if the 
weak parametric downconverter source, described in Ref. 
[ 5j is replaced by a high-flux optical parametric amplifier 
[ 13 . The visibility saturates at 20% in the limit of large 
gain, but this is quite sufficient for some lithography pur- 
poses, as well as for 3D optical holography used for data 
storage. 
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